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HEAT-REJECTION AND WEIGHT CHARACTERISTICS OF FIN- 
TUBE SPACE RADIATORS WITH TAPERED FINS 

By Henry C. Haller, Gordon C. Wesling, 
and Seymour Li el le in 

Levis Research Center 
SUMMARY 

An analysis was made of the radiant blackbody heat-rejection and weight 
characteristics of a direct-condensing radiator with a linearly tapered fin be- 
tween two tubes* Considered in the analysis are such factors as tube armor pro- 
tection, mutual irradiation occurring between the fin and the adjoining tube 
surfaces, and fin base thickness and taper angle. Results of the thermal anal- 
ysis for three fin tapers are presented in terms of a conductance parameter in- 
volving a comparison of the radiation capability of a fin with its heat conduct- 
ing capability, the ratio of fin tip to base thickness, the ratio of tube out- 
side radius to fin lengthy and the ratio of fin base thickness to fin length. 

For practical configurations, the effects of fin base thickness and taper angle 
on total heat transfer are negligible. 

The weight -minimization procedure used in this analysis substantiated the 
fact that the maximum heat rejection per unit weight occurs at conductance pa- 
rameters approximately equal to 1 with tapered as well as with const ant -thickness 
fins. For minimum weight, the tapered fin is thicker at the root and is longer 
than the corresponding minimum- weight constant-thickness fin. The advantages 
gained by tapered fins are highly dependent on the choice of required armor pro- 
tection, material, heat-rejection load, and surface temperature. For the illus- 
trative calculations considered, in which a 1-megawatt- system radiator at 
1700° R with a 0.90 probability of no punctures in 500 days was prescribed, fin 
tapering results in a maximum weight saving of 10 percent for a triangular fin 
as compared to a constant-thickness fin. A similar result was obtained for a 
radiator operating at 1160° R for a 30-kilowatt system. The effect of fin 
tapering on radiator planform area was also investigated. Results indicate that 
an increase in planform area occurs with increasing fin taper at conditions of 
maximum heat rejection per unit weight. These conclusions are based on fin and 
tube heat transfer without consideration of the fluid pressure drop or the 
header weight contributions that may decrease the indicated advantages of taper- 
ing. 


INTRODUCTION 

A major problem in the design of closed-cycle space powerplants is the re- 
jection of large amounts of energy by radiation. Several_ studJLes indicate that 



the meteoroid-protected radiator may he the heaviest component of a space power 
system, and weight optimization is of prime importance for heat-rejection sys- 
tems in space (refs* 1 to 3). The use. of fin surfaces between fluid-carrying 
tubes has been recognized as a practical design concept for reducing radiator 
weight. The characteristics of radiating fin surfaces of several configurations 
have been presented in a number of papers (e. g. , refs. 4 to 6). The investiga- 
tions presented in references 5 and 6 show that linearly tapered fins have 
greater ratios of heat rejection to weight than fins of constant thickness and 
approach the minimum-weight cusp profile analyzed in references 7 and 8. 

The purpose of this report is to analyze the heat-rejection and weight 
characteristics of central-fin-tube radiators incorporating linearly tapered 
fins and to identify potential weight savings over radiators with fins of rec- 
tangular* profile. The analysis considers the radiant interaction between tube 
and fin surfaces and the total weights of both fins and tubes. The first part 
of the analysis applies the principles of reference 9 to heat rejection from a 
fin-tube radiator with linearly tapered central fins. The second part considers 
the overall ratio of heat rejection to radiator weight for several applications 
of interest. The radiator weight includes armor to protect the primary tubes 
from meteoroid penetration as recommended in reference 10, but it does not in- 
clude inlet or outlet headers. 

The heat-rejection analysis is carried out for the situation where the fin 
and the tube are of infinite length in the direction normal to the cross sec- 
tion of the tube. The analysis treats fins having various degrees of linear 
taper, various ratios of tube outside radius to fin half-length, and a range of 
dimensionless conductance parameters. The minimum- weight cusp-profile fin was 
not treated in the analysis presented herein. 
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SYMBOLS 


surface area, sq ft 
vulnerable area, sq ft 
planform area, sq ft 
penetration correction factor 


velocity of sound in material. 



ft/sec 


tube diameter, ft 
Young's modulus, lb/sq ft 

angle factor, fraction of energy leaving surface that is incident 
upon another surface 

acceleration due to gravity, ft/sec 2 
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thermal conductibity,. Btu/(hr) (ft) ( °R) 
half-length of fin, ft 

ctL 2 T-u 

dimensionless conductance parameter, EL 

ky 0 

prohahility of zero punctures 
power plant output, Mw 
heat flow, Btu/hr 
tube radius, ft 
temperature, °R 

average meteoroid velocity, ft/sec 
weight, lb 

normalized distance coordinate, x/L 

coordinate measuring distance along fin, ft 

normalized half-thickness coordinate, y/y Q 

ratio of fin tip to base thickness, y^/y Q 

coordinate measuring half-thickness of fin, ft 

midfin half -thickness, ft 

half-thickness of fin at base, ft 

radiator tube length, ft 

constants in penetration formula 

tube wall thickness, ft 

efficiency or effectiveness 

cycle thermal efficiency 

normalized temperature, t/T^ 

d0 

slope of temperature profile, —L 
approximation to 6^ 
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p 


density, lb/cu ft unless otherwise specified 
a Stefan -Boltzmann constant, 0.1713x10“^ Btu/(hr)(sq ft)(°R 4: ) 

t * exposure time, days 

CD* CD** 

J*jJ**V 3112163 111 fig * 1 

Subscripts: 
a armor 

b tube Base surface 

c tube liner 

cond conduction 

f fin 

i inside 

0 outside 

p particle 

r radiator 

rad radiation 

x coordinate measuring distance along fin, ft 

1 tube surface 1 

2 tube surface 2 

3 adjacent fin surface 


FIN TEMPERATURE PROFILE 
Assumptions 

The analysis considers the general case of a linearly tapered fin of length 
L and taper ratio Y^ attached to a tube with a constant surface temperature 
Tt>. Energy input to the fin comprises heat conduction along the fin from the 
tube -fin interface and incident radiation from the adjacent tapered fin and the 
two tube surfaces. Figure 1 shows the geometrical setup used in the analysis of 
the fin-tube radiator panel, and figure 2 indicates the fin profiles studied. 
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The specific assumptions used in the development of the fin heat-transfer 
relations are 

(1) The radiator surfaces act as hlackhodies with incident and emitted 
radiation governed “by Lambert’s cosine law. 

(2) Incident radiation from external sources is negligible* 

(3) All conditions are constant along the length of the tube* 

(4) Heat is radiated from both upper and lower surfaces of the fin and the 

tube* 

(5) The tube surface is at a constant temperature longitudinally and 
circumferentially. 

( 6) The heat flow in the fins is steady-state one-dimensional flow, and the 
fin base Is at the temperature of the tube* 

(7) Material properties are constant and evaluated at the fin base tempera- 
ture* 


(8) The development of fin and tube angle factors is based on an infinite 
length of tube and fin* 

(9) The fin surface area is based on the length of the fin measured from 
the centerline of the fin at Its base*^- 


Formulation of Equations 


When an element of the fin surface in figure 1 is considered and the pre- 
vious assumptions are employed, the law of energy conservation can be expressed 
as the energy balance between the net heat transfer due to conduction and radia- 
tion: 

^Qcond + ^-Qrad “ 0 (l) 


The net internal heat conduction through the element for half of the fin thick- 
ness is expressed as 


^-Qcond ~ 



( 2 ) 


The net radiant heat rejection from the element of area on the upper surface of 


iThe increase in fin surface area due to an increase in fin thickness Is 
less than one-half of 1 percent of the fin length for a practical limiting case 
given by y 0 /R 0 < 0*1, R^L = 1*0, and R 0 = 1.0 in. 
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the fin is composed of its emission minus the incident energy. The emitted 
energy from the surface element is given hy 


0T 4 

COS 00* 


dx 


(3) 


and the incident energy on the element from the tube surfaces and the opposing 
fin surface is 


/ 2L 

/ T^x 3 -x (4a) 

•'L 

where Aq and Ag are, respectively, the exposed surface area of the tube 
quadrants adjacent to the fin surface in question. 

For practical fin-tube configurations in which y 0 /R 0 is generally not 
greater than 0.10, the amount of energy incident on the fin surface from an op- 
posing fin can "be neglected, since the angle factor from the opposing fin to the 
element on the fin will he considerably less than the angle factors from either 
tube. The net exchange of energy is made additionally insignificant when it is 
noted that the average temperature of the portion of the opposing fin that would 
contribute the most to an element near the fin tip is at the same temperature as 
the opposing fin. When this energy input is neglected and the reciprocity rela- 
tions 


A A-x - z dn# 


' x— 1 


and 


A 2 f 2-x = z 


dx 


COS CD" 


s* r x- 2 


are applied to equation ( 4a ) j the incident energy expression 

4 

OT^Z 

COS CD* (*x-l + 


(4b) 


is produced. The net radiation leaving the surface element is then given by 


4t W = SSV^[ t4 - T ^ (F *-1 + r x-2)] ta (5) 

and the complete energy balance on the element can be written as 

E ( kF f) - P - ’fe-i + f *- 2 >] ( 6) 

This differential equation describing the fin temperature distribution can 
be rewritten in a dimensionless form by introducing the following parameters; 
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Y = 


e = 


JL 



f- = i “ x d " Y z) 

J O 


In -berms of these new variables, equation ( 6) can he written as 


_d_ 

dX 



where 




crL 2 T^ 

k Y 0 


(7a) 


(7b) 


and 


COS CD* = ■ ( 7c) 

V 1 + a - ij ) 2 (^f 

The angle factors in equation (7a) are evaluated by using a relation 
(ref. 11, eq. (31-38)) that applies to parallel surfaces of infinite 
length (i*e. , tuhe and fin length, Z = <») . lor the configuration of figure 1 
the angle factors are 


\ (sin i|r** - sin \|/*) (8a) 

F x _]_ = \ (sin i** - sin \|r*) (8b) 


The evaluation of the angle factors in terms of the parameters X, R q /L, y 0 /L, 
and Yi/Yo is carriel out hy using the geometry and the nomenclature of fig- 
ure 1. The exact expressions developed include the effects of fin thickness and 
fin taper on the angle factors. The required angles that are used in equa- 
tions (8a) and (8h) are 





= sin” 1 


vg? *■)"-»/ 


Rr 


+ tan"! 


+ X 


fy Q 1 - (1 - Yj)x] 

L L x + ^ J 


+ sin 


¥ ^ 

.V 1 + t 1 - (t) 




With the previous formulation for the angle-factor expressions (eqs* (8a) 

and (Ofb))j the exact evaluation of the fin energy balance (eq* (7a)) is complete* 

Calculated variations of the angle factor from adjacent fin to tube, in- 
cluding the effects of fin thickness and taper ratio as obtained from equa- 
tion (8a), are shown in figure 3, As expected, fin taper and thickness reduce 
the angle factor. This comparison indicates that a sizable error in the angle- 
factor from the entire fin to a tube can result when large fin thicknesses are 
not accounted for (as in ref. 9). The true measure of the result of neglecting 
fin thickness and taper angle, however, will depend on the corresponding effects 
on the net heat transfer. 


Computational Procedure 


It was necessary to use numerical techniques to solve equation (7a) for 
and for the G profile from X = 0 to X *= 1* The finite-difference 


^X=Q 

technique used was 
applied to ^ 


a form of Kalaba's method* Central finite differences 


dX 


( T if) 


yield 


± (• 
ax v 


(Ye) , - (Ye) 


Y — ) = 

ax/i 


1 2 


AK 


( 9 ) 


Applying central aifferences to 9 yieias 


(ft 


ax \ ax 4 


1+ l 


i+1 


- ftf) ~ Y ( 0-f - ©i_i) 


i-i 


(AX) : 


( 10 ) 
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The term 0 £ was approximated by a finite Taylor series, namely, 

+ “ ^j[) or of = 40^0? - 30^, w ^ ere % is an approximation to 

6j_+ Thus, 0^ hecomes a linear function of 0^* 

The finite-difference equations were set up in three -diagonal form, were 
reduced to two-diagonal form, and were solved by gaussian backward substitution 
by using an arbitrary initial guess for 0* This process was repeated with the 
new computed values for 0. Four such iterations yielded an accuracy of ±1 in 
the fifth decimal place, or, roughly, each iteration yielded an additional deci- 
mal place. An IBM 7090 electronic digital computer was used for the numerical 
solution. 


Results 

Each solution of the fin energy equation provided a temperature distribu- 
tion along the fin. Results are plotted in figures 4(a), (b), and (c) as func- 
tions of X for several parametric values of N c and Rq/L for triangular, 
trapezoidal, and rectangular fins, respectively. The view factors used to ob- 
tain the curves of figure 4 were based on the simplified case of y 0 /L = 0. 

The effect of the fin thickness parameter y 0 /L on the temperature distribution 
was calculated to be small for the cases shown in figure 4. A maximum error in 
0 of 2 percent results at the midpoint of the fin (X = l) for y 0 /R Q = 0.1, 

Ro/L - 1.0, and large values of N c . For values of N c - 1.0 (near optimum), 
the error in neglecting fin thickness is less than 1 percent at these conditions. 

The curves of figure 4 exhibit the expected trend of small temperature 
drops being associated with low values of N c , and the R 0 /L ratio contributes 
only a minor effect at this condition. As U c increases, so does the fin tem- 
perature drop and the importance of the ratio R 0 /L. A comparison of figures 
4(a), (b), and (c) shows that increasing the fin taper increases the temperature 
drop along the fin for fixed values of the conductance parameter and the ratio 
R q /L. This is plausible since the fin cross-sectional area at any point on the 
fin is decreased as the taper is increased, and thus the thermal resistance is 
increased. 


RADIATOR HEAT TRANSFER 


Fin Heat Rejection 


After the temperature distribution and the slope of the temperature dis- 
tribution curve at X = 0 have been determined, the net heat transferred by the 
fin can be calculated. The net heat loss from one side of the fin of length L 
is equivalent to the amount of heat conducted into half of the fin base thickness 
(fig. 1). Fourier’s he at- conduction equation is 


Of = - ky Q Z 



- kz r 0 T -b /ae\ 

L \3X/x=0 


( 11 ) 
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This equation can he placed in dimensionless form and can he evaluated by using 
the solutions of the fin energy equation (7a). The comparison of this fin 
energy rejection to the total heat loss from one side of a "black isothermal fin 
of length L can he expressed as 



Equation (12) is the definition for the efficiency of the fin. 

For the case of N c = 0 (infinite thermal conductivity and constant tem- 
perature), equation (12) is indeterminate. For these conditions, however, the 
fin efficiency can he found from simple geometrical considerations. Since the 
heat loss is the difference between emitted and incident energy, the fin effi- 
ciency becomes 


Qf* 

OLZT^ 


= 1 



( F X-1 + F X-2) dX 


(13) 


Thus, when R 0 /L = "Hf = !• and when R 0 /L *+ », % “*• 0. 

Results of the fin efficiency calculations are shown in figure 5 for varia- 
tions in N c , Rq/L, and fin taper for the simplifying case of y 0 /L = 0. In- 
creasing the conductance parameter N c or the fin taper y^/y Q decreases the 
fin efficiency because of the corresponding reductions in the fin temperature 
profile (fig. 4). For the limiting case when R 0 /L = 0, there is no radiant 
interchange between tube and fin. This is the plane fin case that is treated in 
references 4 and 6. The decrease in rjf with increasing R 0 /l> is a reflection 
of the increasing effects of tube occlusion. 


If the fin thickness ratio is considered in the fin heat rejection, the fin 
efficiency as described by equation (12) will increase. This is caused by a re- 
duction in the occlusion of the tube on the fin, which permits a larger amount 
of energy to escape. Numerical results indicate less than a 2-percent error is 
incurred in T) f when y 0 /L is neglected for the case of R q /L = 1.0, 
y 0 /Ro = 0.1, and N c in the range of interest. 


Tube Heat Rejection 

The net heat loss from the tube surface is equal to the energy emitted by 
the exposed tube surface minus the incident radiant flux, which is composed of 
emission from the fins and the opposing tube surfaces. If one-quarter of the 
surface of tube 1 in figure 1 is considered, the tube-surface heat loss per unit 
length of tube can be written as 
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Inasmuch as the incident energy from the fin, at X = 0 to X = 1, to tube 2 is 
equal to the incident energy from the fin, at X = 1 to X = 2, to tube 1, the 
fin-to-tube angle factor in equation (14) can be expressed as 




(*x-l + F x _ 2 )dx 


(15) 


The determination of the angle factor from tube 2 to tube 1, including the ef- 
fect of fin thickness, is obtained from angle-factor algebra and is written as 


F 2 _i = 1 - 


I - sin" 1 


(Vo) 

w 


1 + 


cos as 




— f 


( F x-1 + F x _ 2 )dx 


(16) 


Substituting eqLuations (15) and (16) into equation (14), nondimens ionali zing, 
and comparing the actual heat loss % with the ideal emission from the 
quarter-circle surface of the tube yield 


% 


Z 




COS CD* COS CD* 


f 

*r\ 


(Fx-l + f x _ 2 )(2 - 04) ax 


•= % 


(17) 


This expression, which is the tube-heat-rejection efficiency, is identical in 
form to the equation presented in reference 9 with the exception of the cos cd* 
term, which is a function of the fin taper (eq. (7c)). The effects of fin taper 
and fin thickness are also introduced in the angle-factor terms and F x _ 2 

obtained from equations (8a) and (b) and from the numerical values of the fin 
temperature profile e(x) obtained from the solutions of equation (7a). 

When the conductance parameter approaches zero, the fin is isothermal, and 
equation (17) reduces to 
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The tube efficiency for this case, however, is not affected hy fin taper. The 
cos 03* term and the angle-factor expressions, which contain the fin taper and 
thickness parameters, cancel out numerically for y 0 /R 0 equal to or less 
than 0. 1. 

Figure 6 shows a plot of the tube efficiency as a function of the ratio 
R 0 /L for selected values of N c and fin taper when the effects of fin thick- 
ness, y 0 /L, are neglected. For the N c = 0 curve, the R 0 /L = 0 condition 
represents the case of a tube with an infinite isothermal central fin. For this 
situation, the fin occlusion reduces the tube efficiency to 0.818. The other 
limiting condition of R 0 /L oo represents the case of isothermal touching 
tubes, for which the efficiency approaches 2/jt. Tube efficiency increases as 
the conductance parameter of the fin gets larger, since the fin emits less to 
the tube because of its reduced temperatures. Fin taper also increases the tube 
efficiency because tapering reduces the fin temperature profile. As discussed 
in reference 9, however, these curves should approach the single limiting values 
of the curve for N c = 0 at R 0 /L = 0 and oo. The dashed lines at low values 
of R 0 /L in figure 6 are faired curves. 

Consideration of the effects of fin thickness on the tube efficiency indi- 
cates that a reduction in efficiency occurs when the fin thickness is increased. 
This is caused by a reduction in the prime surface area of the tube. The reduc- 
tion in r)^ is less than 2 percent for the case when the ratio of fin thickness 
to length is 0.1, R 0 /L = 1* and Nc = 1* Even though the tube area reduction 
is approximately 10 percent at this condition, the effect on the tube efficiency 
is small since the tube area near the fin contributes little net emission to 
space because of the occluding effect of the fin. 


Total Fin-Tube Heat Rejection 

It is desirable to compare the tapered-fin and rectangular-fin radiators on 
a common basis, that is, the same R 0 /L and conductance parameter. It is also 
important that a fin-tube section of both radiators be compared since variations 
in fin length and base thickness have marked effects on the heat rejection from 
the tube or the fin when treated separately. If the definitions of fin effi- 
ciency in equation (12) and of tube efficiency in equation (17) are used, the 
following useful definition of overall fin-tube effectiveness is formulated: 


1r = ^4 + £) 

Equation (19) is for a quarter section of a tube and a fin, but it is identical 
to that for an entire fin-tube section of length 2L(l + R 0 /L) radiating from 
two sides. 

Figure 7 shows the effectiveness of the entire fin-tube radiator that was 
obtained by using equation (19) for the case of y 0 /L = 0. The trends with the 
taper parameter Y i and the conductance parameter N c are in the same direc- 


^f + % 


l + 2o 

L 


it R o 


2 L 


(19) 
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tion as those exhibited by the fin efficiency curves* Factors that reduce the 
fin thermal resistance improve the overall fin-tube effectiveness. The greatest 
improvement in radiating effectiveness occurs with increasing R 0 /L because 
more of the radiator is at the base temperature. 

The total effectiveness of the fin-tube radiator is equal to 1 for the two 
limiting cases of R 0 /L = oo and the conductance parameter N c = 0. In both of 
these cases the entire fin and tube section is isothermal and acts as a black- 
body radiator. The total effectiveness approaches the fin efficiency for the 
third limiting case, R 0 /li = 0. At this condition, the tube does not contribute 
to the total radiator effectiveness even though the tube efficiency is finite, 
as shown in figure 6. 

The inclusion of the fin thickness parameter y 0 /L resulted in less than a 
1-percent variation in total fin-tube effectiveness for y Q / R 0 =0.1 at 
R 0 /L = 1 and conductance parameters of interest. Thus, within the range of 
conditions calculated, neglecting the effects of fin thickness on the total 
radiator heat transfer can generally be justified. 

It is also of interest to know what fraction of the overall heat loss is 
transferred from the fin. A plot of the ratio of fin heat rejection to total 
radiator heat rejection is presented in figure 8. Inspection of the figure re- 
veals that the fraction of heat lost by the fin is greatest for small values of 
R 0 /b and N c . The tendency of Qf/Q r “to increase with decreasing R 0 /L is 
related to the increase in the fin surface area relative to the tube surface 
area. On the other hand, the changes in Qf/Qr with N c and fin taper are due 
to changes in the fin temperature distribution. 


RADIATOR WEIGHT 

The heat-rejection analysis of the fin and the tube was nondimens ionalized 
so that the results could be used for general design purposes. Using the pre- 
vious results for making a decision as to the merits of employing tapered fins 
in radiator designs is impractical without consideration of the total weight of 
the fins and the tubes. It is necessary, therefore, to consider the ratio of 
heat rejection per unit weight Q/W. 

Input information required from the heat-rejection analysis consists of the 
overall radiator effectiveness (eq. (19)) as a function of Nc^ Ro/Lj Yi/y q > 
and y 0 /L. Because of the additional variables of tube internal diameter, 
meteoroid shield thickness, radiator materials, radiator temperature level, and 
system power level, weight and area optimizations can be made only for specific 
cases. 

For comparison purposes, it was considered sufficient for the panel geom- 
etry to match only the total required heat flow. Consideration of fluid pres- 
sure drop and vapor and liquid headers was not included in the formulation of 
the radiator weight relations. This simplification in essence permits the 
analysis of a single tube and fin without regard to the number of tubes or to 
the subsequent individual tube length. A tube-to-tube span is considered in 
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the thermal analysis, however, so that the radiation interchange between fin 
and tube is properly included. 


Weight Ratio 


The heat rejection per unit weight of a fin-tube radiator is given by 


£ = Q/Z 

w w/z 


( 20 ) 


where, for radiation from both surfaces, 

I ' ^ + £)%< < 21 > 

For a radiator tube composed of a thin inner liner of thickness 5 C and an 
armor sleeve of thickness 5 a , the outside tube diameter may be expressed as 

Do = D i + 2S C + 26 a (22) 


The fin and tube weight per unit length of tube can be calculated by sum- 
ming the individual weights of the fin, the tube liner, and the tube armor to 
obtain 


£ « Pcfpi 


(23) 


2 3 

When y 0 = cLTg/kN c and D 0 from equation (22) are introduced into equa- 
tion (23), the equation for weight per unit length becomes 


W 

Z 


p c Jt5 c (Di + 5 C ) + P a JcS a (D 1 + 26 c 


+ 8 a ) + 


2p f dT^(l + Y z ) 
kW c 




(24) 

Solution of equation (24) for fixed values of HT C , Tj-,, and Y i for a given 
radiator requires input values for the liner and armor thickness. The former is 
generally prescribed as some function of the tube inside diameter, while the 
latter depends on the meteoroid protection criterion selected. The tube liner 
thickness is arbitrarily prescribed as 

6 C = 0.04 Di (25) 

The meteoroid protection criterion used in this analysis is that of refer- 
ence 10 with revised inputs based on reference 12, as recommended by the authors 
of reference 10 in a personal communication. According to reference 10, the re- 
sultant equation for the armor thickness is 
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where (3 = 1.34, a = 1.75, Vp = 98,400 feet per second, pp = 0.44 gram per cubic 
centimeter, and a. = 0. 5 3x10“ 10 g 3 P/(sq. ft) (day) (Whipple value without Earth 
shielding) . The radiator vulnerable area, taken as the exposed surface area of 
the tube, is given "by 


or, from equation (21), 



(27) 


(28) 


The heat-rejection load is related to the cycle power output and overall effi- 
ciency by 

Q r = 3. 41 3x10 6 & e (~~ - l) (29) 

Thus, an iteration solution for w/z involving equations (19), (26), (28)^ 
and (29) is required. 


Calculations 

Two specific cases are given in order to show the effect of power level, 
temperature level, and material specification on the radiator heat rejection per 
unit weight. The first case is a 1-megawatt electrical output powerplant with a 
15-percent cycle efficiency and a 1700° R radiator temperature. The tube armor 
and the fin were assumed to be beryllium, and the tube liner was assumed to be 
columbium. The second case is a 30-kilowatt system radiating at 1160° R with a 
15-percent cycle efficiency. For this case the tube armor and the fin were taken 
to be aluminum, and the tube liner was taken to be stainless steel. Radiator 
tube inside diameters of l/2, 3/4, and 1 inch were used for both cases investi- 
gated. A 500-day mission time and a probability of no puncture P(0) of 0.90 
were chosen for the calculation of meteoroid protection thickness. 

Figures 9 and 10 illustrate the variation of the ratio of total heat re- 
jected to radiator panel weight with fin-tube ratio L/Ro and the conductance 
parameter N c for the two examples. The effect of fin tapering in increasing 
q/W is most pronounced for large values of L/Ro, where the fin contribution to 
the total heat transfer is greatest. A peak value of Q/W is attained in all 
cases. Similar curves are obtained for tube inside diameters of l/2 and 1 inch. 

Plots of maximum values of Q/W against L/R 0 for the values of iT c con- 
sidered are shown in figures 11 and 12 for all three tube inside diameters and 
fin tapers. At both power levels, fin tapering results in an increase in peak 
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q/W of approximately 10 percent for the triangular profile (Yj = 0). From 
studies of fins alone (ref. 7), gains of up to 30 percent are indicated. The 
difference in the case of the radiator panel is explained hy noting- that the 
armored radiator tubes account for a sizable portion of the total -weight and heat 
rejection. Thus, the effect of fin tapering on total weight is effectively re- 
duced for the fin and tube configuration. 

Peak q/w occurs at a value of N c of approximately 1 in all cases con- 
sidered. The corresponding values of L/r q for peak q/w vary somewhat with 
fin taper and tube diameter at a given power and temperature level. Fin taper- 
ing tends to increase the value of L/R 0 for peak Q/V. For the 1-megawatt 
system, the best L/R q for maximum q/W varies from approximately 2.5 to 4.0. 
For the 30-kilowatt system, the best L/R q increases to values from approxi- 
mately 5.0 to 8.5. 

The 10-percent weight saving indicated in the two examples must be quali- 
fied, in that the effects of pressure drop and vapor and liquid header weights 
were not considered. In effect, the optimum triangular fin is longer than the 
optimum rectangular fin. Thus, some of the potential weight savings might be 
offset because the headers for the triangular-fin radiator would be longer and 
heavier. In this comparison the number of tubes and their length would have to 
be assumed constant and would be obtained from pressure-drop considerations in 
the tubes and the headers. 

Consideration of the fin-thickness parameter y 0 /L in the weight- 
optimization procedure proved to have very little effect on the heat rejected 
per unit weight in the vicinity of the design point (N c = 1.0). At some ex- 
treme nonoptimum points, a maximum difference of 3 percent or less was obtained. 


RADIATOR GEOMETRY 


In addition to radiator fin and tube weight, the required planform area of 
the radiator panel is, in many cases, also of considerable interest with respect 
to proper integration of the radiator and the space vehicle. Radiator planform 
area is obtained from equation (21) as 


A* = Z(D 0 + 2L) 


Qr 


(30) 


Thus, planform area will vary inversely with total fin-tube effectiveness, or, 
from figure 7, the area will generally increase with increasing L/Ro. 

Figures 13 and 14 illustrate the calculated variations of planform area 
with L/r 0 , Re* and fin taper for the two power cycles for peak Q/W conditions. 
Fin tapering increases the required planform area for a given value of N c . 
Similar results are obtained for the other inside tube diameters. 

Other factors of interest with respect to the geometry of the radiator are 
the fin base thickness and the tube armor thickness. The effect of fin tapering 
on these factors for a l/2-inch tube inside diameter is illustrated in figures 
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15 and 16 for peak Q/W conditions. In 'both cases, fin tapering results in 
sizable increases in base thickness and minor reductions in tube armor thickness. 
Increased fin base thickness may be of advantage with respect to joining or 
fabrication considerations. 


SUMMARY OF RESULTS 

From an analysis of the heat-transfer characteristics of fin-tube space 
radiators with linearly tapered fins (rectangular, trapezoidal, and triangular), 
it was determined that 

1* Fin-tube radiators with tapered fins have lower thermal efficiencies and 
fin temperature profiles than radiators with rectangular fins. 

2. Although fin base thickness and taper angle can have a marked effect on 
the angle factors between fin and tube, their net effect on the total heat re- 
jection is relatively minor. 

Sample calculations of tapered-fin radiators for typical 1-megawatt and 
30-kilowatt power systems showed that 

1. Tapered fin radiators can be about 10 percent lighter than radiators 
with fins of constant thickness for the assumptions used for the power cycle and 
meteoroid protection. The calculation, however, did not include the effects of 
header weights and fluid pressure drop. 

2. Conductance parameters of approximately unity resulted in minimum fin- 
tube weight for all tapers considered. Corresponding values of the ratio of fin 
length to tube outside radius varied from 2.5 to 4.0 for the 1-megawatt system 
to 5.0 to 8.5 for the 30-kilowatt system. 

3. Tapered fins are longer and have greater base thicknesses than fins of 
constant thickness for maximum heat rejection per unit weight. 

4. Fin tapering tends to increase radiator planform area at maximum heat 
rejection per unit weight. Planform area, however, varies most markedly with 
the ratio of fin length to tube outside radius. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, October 29, 1963 
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(a) Rectangular; ratio of fin tip to base 
thickness, 1- 0. 



(b) Trapezoidal; ratio of fin tip to base 
thickness, 0. 5. 



(c) Triangular; ratio of fin tip to base 
thickness, 0. 

Figure 2. - Fin profiles studied. 


Angle factor, Ff_i 










Normalized temperature 
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Normalized distance, X 

(c) Rectangular fin; fin taper ratio, y^/y 0 > 1.0. 

Figure 4. - Concluded. Fin temperature distribution. Fin thickness ratio, y 0 /L, 0. 


Pin efficiency, = ^/aT^LZ 



Figure 5. - Fin efficiency as function of ratio of tube radius to fin half-length. 
Fin thickness ratio , y Q /L, 0. 
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Ratio of tube radius to fin half-length, R 0 /L 

Figure 6. - Tube efficiency as function of ratio of tube radius to fin half-length. 
Fin thickness ratio, y 0 /L, 0, 






Figure 7. - Total fin-tube effectiveness as function of ratio of tube radius to fin 
half-length. Fin thickness ratio, y 0 /L, 0. 


27 




Heat rejection per unit weight, Q/W, Btu/(hr)(lb) 


2 * 

10X10 3 

8 


4 


2 

8xl0 3 
6 ■ 



4 


1 


Pigu: 

poi 

ari 





Heat rejection per unit weight, Q/w, Btu/(hr) (lh) 





Maximum heat rejection per unit weight, Q/W, Btu/(hr) (lb) 




Maximum heat rejection per unit weight, Q/W, Btu/(hr)(lb) 





Planform area 
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Figure 13. - Radiator planform area requirements for 1-megawatt 
system radiating at 1700° R. Beryllium fins and tube armor; 
tube inside diameter, l/2 inch. 




Fin base thickness, 2y 0 , ft Armor thickness 




Fin base thickness, 2y Q , ft Armor thickness 




Figure 16. - Radiator physical dimensions for 30-kilowatt power system 
radiating at 1160° R. Aluminum fins and tube armor; tube inside 
diameter, l/2 inch. 
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